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Introduction 

In this work, we are concerned with the well studied advection-diftusion equation in two 
space dimensions which we write in the following form: 



where u is a given incompressible velocity field, k is a given diffusion constant, c is the 
unknown concentration function and the factor of two is included for future algebraic 
convenience. Most of the difficulties encountered in the study of this equation arise from 
the fact that the advection operator 2u- V and the diffusion operator kA exhibit different 
symmetries: The advection operator is skew-symmetric while the diffusion operator is 
symmetric and hence self-adjoint. As a result, the combined advection-diffusion operator 
possesses, in general, neither symmetry. This lack of symmetry complicates the study 
of solutions of Eq. (U especially in the dynamically interesting case of small diffusivity 



For a very particular class of advection-diffusion operators, namely those associated 
with irrotational (potential) velocity fields u, it is well-known that a simple point 
transform maps ([I]) into a selfadjoint Schrodinger-like problem [2]. Consequently, 
the transformed equation can be studied using standard techniques from quantum 
mechanics, in particular WKB approximations, and the results have important 
implications for the original equation. This approach has been used in the past in 
the context of Fokker-Planck equations [3]. From both the physical and dynamical- 
systems points of view, cases where the velocity field is time-dependent are of particular 
interest. In one spatial dimension, all flows are potential and this fact allows one to 
apply point transforms to fairly general one-dimensional non-autonomous convection- 
dominated parabolic equations. In doing so, difficulties arising from the fact that 
convection-dominated parabolic equations do not satisfy the spectral-gap condition, 
required for all currently available proofs of the existence of inertial manifolds, are 
completely circumvented. The existence of inertial manifolds, in turn, allows for a 
complete decomposition of the convection-diffusion operator in terms of Floquet bundles. 
Using averaging techniques and inverse scattering theory, it is possible to further 
transform the non-autonomous parabolic equation, at least in one spatial dimension, 
into an autonomous Schrodinger-like equation [3]. 

Point transformations of this kind prove useful not only for linear autonomous and 
non-autonomous parabolic equations, but also for certain nonlinear parabolic equations. 
Point transforms have been used to circumvent the spectral-gap condition and prove the 
existence of inertial manifolds for a class of nonlinear nonlocal Fokker-Planck equations 
and for a class of viscous Burgers equations with low-wavenumber instability in both 
one and two space dimensions [5, 6J. 

For the latter, the necessary point transformation is the well-known Cole-Hopf 
transformation. The existence of inertial manifolds implies, in particular, that possibly 
very complicated global attractors can be embedded in smooth finite-dimensional 




(1) 
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manifolds on which the original PDE reduces to a finite system of ODEs. As such, 
the finite-dimensionality of the global dynamics is rigorously established. 

These remarkable results for equations featuring advection by potential flows 
naturally lead to a search for an extension of the point transform to equations featuring 
advection by rotational flows. As we show here, no exists no simple point transform 
that maps Eq. [1] into a Schrodinger-like equation. In what follows, we propose the 
combination of a point transform with an appropriate Lie-transform in order to obtain 
a similar result. The approach is perturbative in the diffusivity and leads to a rather 
intricate system of nonlinear 'balance' equations that the sought-after transformation 
must satisfy. A complete understanding of this system requires methods of nonlinear 
functional analysis and this continuing work goes beyond the scope of the present 
paper. Our purpose here is simply to provide the kernel of the combined point-operator 
transform and illustrate its applicability in the context of a simple, canonical example. 



Point and Lie Transforms 



We begin this section by assuming that u is a given, not necessarily incompressible flow, 
and illustrate how one arrives at the point transform described in the introduction for 
irrotational flows, and why the same approach fails for flows with non-vanishing curl. 
For simplicity, let us assume that the velocity field does not depend on time while noting 
that all results can be extended to the time- dependent case. Consider a simple point 
transform of the form 

c = e^ /K v . (2) 

Substitution of ((2]) into ([I]) yields the following equation for v: 

v t + Bv = —fv + kAv . (3) 

K 

The operator (1/ n)f + kA on the right hand side of the transformed equation is again 
symmetric, and the potential / given by 

/ = <\> 2 X + 0j - 2u X (f> x - 2u 2 <\> y = | V0| 2 - 2u ■ V0 (4) 

The new advection operator B is given by 

B = 2( Ul - <j) x )d x + 2{u 2 - <f> y )dy - A0 , (5) 

and one can easily verify that B is also skew- symmetric. Now, we can see directly 
that, if u is a potential flow, we can choose <fi such that V0 = u, and we obtain B = 0. 
Thus, we have transformed the advection-diffusion equation involving a skew-symmetric 
advection operator on the left hand side and a symmetric diffusion operator on the right 
hand side into a self-adjoint Schrodinger-type problem given by 

v t = K{- 2 f + A)v. (6) 
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The above calculation shows also that, for a velocity field with non- vanishing vorticity, 
a simple point transform (jSJ) cannot be sufficient. The fundamentally new idea of this 
work is to use a further transform, more precisely a Lie transform, 

v = e K V L + = L (7) 

and to choose the function <p and the operator L in a way that, at least for small 
diffusivities k, the resulting problem becomes self-adjoint. One motivation for choosing 
a Lie transform is the fact that this allows us to systematically construct approximate 
solutions in the limit k — > as we can use the Baker-Campbell-Hausdorff formula [7j to 
find the following expansion up to the order 0(k 2 ): 

w t + Bw + k[B, L}w = ( -/ + [f, L] + | [[/, L], L^j w + kAw (8) 

Note that, at order 0(k 2 ) and higher, a variety of terms appears, including commutators 
that involve not only /, B and L, but also A and L. Our aim is to choose <fi and L in 
a manner that the equation for w involves only self-adjoint operators. Since B is skew- 
symmetric and L is symmetric, one can easily verify that [B, L] is symmetric. On the 
other hand, [/, L] will be skew-symmetric and [[/, L], L] will be symmetric. Therefore, at 
leading order in the expansion using k as small parameter, in order to balance the skew- 
symmetric operators in a manner that their influence on the evolution of w vanishes, 
we require 

B = [f, L] . (9) 

Assuming that we were able to find a and an operator L such that the balance ([9]) of 
gradients is satisfied, the resulting equation for w becomes 

w t = K Aw, A=-hB,L]+A + \f, (10) 

where we have omitted higher order terms. Note that, in general, the commutator [B, L] 
will not vanish and A will not be a Schrodinger operator. By construction, however, 
[B, L] is symmetric and hence A will be symmetric. In order to satisfy symmetry 
requirements at higher orders in the expansion in k, it is also possible to expand and 
L in an asymptotic series as 

<f> = Yl L = Y1 KkLk 

k k 

in order to derive a hierarchy of equations for 0& and In this paper, however, we 
will only consider expansions up to 0(k 2 ). 

So far, the only requirement we have imposed on the operator L is its symmetry, 
L = L + . Note that, at this stage, we have not imposed any conditions on the function 
<j), so that there is great freedom in the choice of L and (f). The assumption L = L + can 
also be modified, to e.g. L = —L + , or dropped entirely by writing L = L s + L a where L s 
denotes the symmetric and L a the skew-symmetric part. This will change the condition 
for balancing the gradients in (Q. In this paper, we are restricting our attention to the 
symmetric case. 
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As B is a first-order differential operator and / is a function, it is reasonable 
to search for an appropriate L in the family of symmetric second-order differential 
operators. Symmetry requires L to be of the form 

L = c u d xx + 2c 12 d xy + c 22 d yy + hd x + b 2 d y (11) 

with the conditions 

b\ = c nx + c 12y , b 2 = c 12x + c 22y (12) 

One major advantage of choosing a differential operator is that the commutators are 
relatively easy to calculate. From (|9]), the balance equations become 

cnf x + c 12 f y = 4>x - Ui (13) 

c\ 2 f y + c 22 f y = <p y -u 2 . (14) 
Strictly speaking, from (Q we actually obtain a third equation given by 

Lf = A<f>. (15) 

An easy calculation, however, shows that if ( IT31T41) are satisfied, ( TT5l) will be satisfied 
as well, and the search for the desired transformation reduces to solving the balance 
equations (1131)141) in conjunction with the equation (j4j) for calculating the potential / 
from the potential <p. Note that solving fll3|14(4l) for a given velocity field u means 
finding functions <f), Cn, Ci 2 , c 22 . The balance equations might appear simple at first 
glance; however they are quite difficult since they are nonlinear due to the fact that / 
depends on through (111) . For purposes of this paper, we shall content ourself with 
providing an example in which the system has a solution which is easy to determine. 

Simple shear flow 

In this section, we will examine what the equations of the previous section mean for the 
simplest example of advection by an irrotational flow: simple shear flow. To summarize 
the previous section, we have shown that, if a solution of the balance equations ( JT3lfT4l) 
together with (@| can be found, the original advection-diffusion equation can be 
transformed (up to terms of order 0(k 2 )) into a selfadjoint problem given by ( JTUI) . Let 
us consider a velocity field of the form 



«= Z (16) 




We assume a ^ (3, so that V x m ^ is satisfied. As u is formed by terms that are 
linear in x and y, we assume a quadratic form for 0: 

<p = 2~ ax2 + CX V + 2~by 2 (17) 
with constants a, b, c to be determined later. Then the expression for / yields 

/ = ( a 2 + c 2 - 2f3c)x 2 + (b 2 + c 2 - 2ac)y 2 + 2{ac + be - era - (5b)xy . (18) 
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Due to the particular structure of the velocity field, we can write 

V/ = f(*V F= ( F " F A ( 19 ) 



with a symmetric matrix F whose coefficient can be directly read from the explicit 
formula for /. The system of equations ( lll3][T4l can now be written in matrix form as 

c *="==( c vr) 

and thus we can express C in terms of {a,b,c} together with {a, (3}. By solving this 
system explicitly, we find that the requirement a ^ (3, yields the condition 

c 2 = ab. (21) 

The explicit formula for C under this assumption is 

2 1q I a2<x + a /3b — /3cb — abc (a + /3)c 2 — 2a/3c \ 
' (a + (3)c 2 - 2a(3c (3 2 b + a/3a - aca - (3ac ' 

where we used the abbreviation 7 = aa — fib. For the explicit calculation of A governing 
the evolution of w, we find that the coefficients of the commutator [B,L]/2 are given 
by the matrix 2MC. Writing 

A - ^ [B, L] = w u d xx + 2w 12 d x d y + w 22 dyy (23) 

we can easily find all the entries of the corresponding matrix W. Direct calculation 
shows that the determinant of W is given by 

a 2 8 2 

det(W) = J- (24) 

hence the resulting operator W is either hyperbolic or parabolic, but not elliptic. 

If one of the constants, a or f3 vanishes, the velocity field reduces to a simple shear 
flow. Without loss of generality, let us consider 

u=(ay,0). (25) 

For this case, we find the matrix as 

1 cb b c 

11 2a 2cra 2 ' 12 2cra' " 2aa 

and 

[B,L] = 2(\--) d xx + 4-d xy . (27) 

Remember that the evolution of w is governed by the operator A defined in (flOj) . For 
the first two terms in the definition of A we have found 

- 1 [B, L] + A = -d xx + d yy - 2-d x d y (28) 

In this case, it follows from (J24l) that the evolution of w is governed by a parabolic 
equation. The potential / is given by 

/ = ( a 2 + c 2 )x 2 + 2{ac + cb - aa)xy + (c 2 + b 2 - 2ac)y 2 (29) 
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Comparison of explicit solutions 

It is instructive to see in detail how the transformation works for a particular initial 
condition, for which both the original equation (0Q) as well as the transformed equation 
can be solved explicitly via the Fourier transform on a finite time interval. This will 
allow us to compare the two solutions and verify the validity of our method. Let us look 
for the solution to the original advection-diffusion equation ([1]) for the shear flow with 
an initial condition given by 

c(t = 0,x,y) = 5(x)6(y). (30) 

Note that, in this case, the analytical solution for c(t,x,y) is well-known and can be 
obtained directly from (JT|) using Fourier transform and the method of characteristics. 
In order to obtain the same solution using the sequence of the proposed transforms, let 
us first set (a,b,c) = (a, 0,0), so that the condition fl2T|) is trivially satisfied. In this 
case, we obtain for the function <fi and the operator L simply 

*=l ax2 > L = i& < 31) 

The transformation of the initial conditions is trivial as 

w(0, x, y) = e- KL v(0, x, y) = e~ KL e-' l,/K 5(x)5(y) (32) 

In order to make this equation meaningful, however, it seems reasonable to assume 
a < 0. In this way, — kL will be regularizing. In frequency domain, we obtain then for 
w(0,u,rj) obviously 

w(0,u,ri) =e K " 2/(2a) (33) 

which decays, if a < 0. The next step consists in solving the equation for w that 
simplifies to 

w t = (^dyy H — (a 2 x 2 — 2aaxy) S j w (34) 

Using a Fourier transform only in y and characteristics, we find the explicit solution of 
this equation as 

W(t, X, 7]) = e ~^ 2 tWxH/K+2 aal x^ + Ao?a^^/{ZK)^^ % ^ _ 2a ixt / k) (35) 

Setting (3 = —k/(2o) > 0, we can write the initial condition for w as 

w{0,x, V ) = -±=e« /, 'M (36) 
z\/7rp 

and obtain the explicit formula for the Fourier transform of the solution, w = w(t, u, rj), 



as 



Mt u rj) = —L= e -("-2aa^)7(47) e -^ 2 * (37) 



with the abbreviation 



7 = Y- ( 1 + 2at ( 1 + JaV ) ) . (38) 
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Note here that the Fourier transform method we employ here is valid only as long as 
7 > 0, i.e., on a finite time interval. After that, the Fourier transform ceases to exist, 
and the method leading to the explicit solution breaks down. We can interpret this as a 
loss of regularity in the frequency domain in this particular step of the method, which 
will actually be gained back in the subsequent step of the transformation. Moreover, 
7 > is only possible given our initial assumption of a < 0. We now transform back 
from w to v in Fourier space via 

v(t, u, rj) = w(t, to, v )e- KU>2/ ( 2a) (39) 

In order to obtain the solution c(t,x,y), we only need to transform the solution (1391) 
back to (x, j/)-space and multiply by exp(ax 2 /(2/t). It is precisely in this multiplication, 
where we get back the necessary regularity (remember a < 0). Moreover, using the 
explicit formulas for Gaussian integrals, one can easily check that, at the end of the 
calculation, the parameter a drops out, as it should since a was only introduced in the 
transformation and not present in the original advection-diffusion equation. 

For the simple example of spatially homogeneous vorticity, L and [B, L] are 
differential operators with constant coefficients and all possible higher order terms in 
([8]) vanish. Therefore, in this particular case, the transformation is exact and holds for 
any value of k. 
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